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ABSTRACT 

We study the damping of the gravitational radiation-driven /-mode instability in ro- 
tating neutron stars by nonlinear bulk viscosity in the so-called supra-thermal regime. 
In this regime the dissipative action of bulk viscosity is known to be enhanced as a re- 
sult of nonlinear contributions with respect to the oscillation amplitude. Our analysis 
of the /-mode instability is based on a time-domain code that evolves linear perturba- 
tions of rapidly rotating polytropic neutron star models. The extracted mode frequency 
and eigenfunctions are subsequently used in standard energy integrals for the gravita- 
tional wave growth and viscous damping. We find that nonlinear bulk viscosity has a 
moderate impact on the size of the /-mode instability window, becoming an important 
factor and saturating the mode's growth at a relatively large oscillation amplitude. 
We show similarly that nonlinear bulk viscosity leads to a rather high saturation am- 
plitude even for the r-mode instability. In addition, we show that the action of bulk 
viscosity can be significantly mitigated by the presence of superfluidity in neutron star 
matter. Apart from revising the /-mode's instability window we provide results on the 
mode's gravitational wave detectability. Considering an /-mode-unstable neutron star 
located in the Virgo cluster and assuming a mode amplitude at the level allowed by 
bulk viscosity, we find that the emitted gravitational wave signal could be detectable 
by advanced ground-based detectors such as Advanced LIGO /Virgo and the Einstein 
Telescope. 
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1 INTRODUCTION 

Neutron stars are exciting cosmic laboratories of matter at 
extreme conditions. However, unveiling the physics of these 
ultra-dense and relativistic systems is a real challenge even 
after more than four decades of astronomical observations. 
The situation is likely to improve significantly with the 
advent of gravitational wave astronomy. The forthcoming 
generation of improved gravitational-wave detectors, such 
as the Advanced LIGO and Virgo interferometers and the 
planned Einstein Telescope (ET), should be sensitive to var- 
ious aspects of neutron star dynamics. The observation of 
neutron star oscillations (normal modes), usually dubbed 
neutron star "asteroseismology" , could become an excel- 
lent t ool for probing the interior properties of these ob- 
jects (|Andersson fc Kokkota3ll998l '). 

Of particular importance for neutron star asteroseis- 
mology are the oscillation modes that can become unsta- 
ble (i.e. grow in amplitude) under the emission of gravi- 
tational radiation. This is the well-known CFS instability 



which develops when a retrogade mode (with respect to the 
stellar rotation) becomes prograde, as a result of rotational 
draggi ng, and then grows while emitting gravitational radi- 
ation (|ChandragekhaJll970l : iRiedman fc Schut j|l978t ). Ex- 
haustive work on the subject has shown that among the var- 
ious stellar modes the ones that could become CFS-unstable 
in realistic neutrons stars, and likely to be strong sources of 
gravitational waves, are the inertial r -mode and the funda- 
mental /-mode fsee lAnderssonI (|2003t ) for a review). 

The study of dissipative mechanisms that could limit 
or entirely suppress the gravitational wave instability of the 
r- and /-modes h as become a small industry si nce the early 
days of this field (|Andersson fc KokkotadboOll ). Given that 
astrophysical neutron stars are viscous systems, a clear un- 
derstanding of the various dissipative mechanisms is a nec- 
essary ingredient of any realistic theoretical model. In fact, 
this is another facet of neutron star asteroseismology as the 
absence of unstable modes can be related to the strength of 
dissipation, thus providing indirect information about the 
stellar interior. For instance, this possibility has been re- 
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cently illustrated by Andersson fc Haskelll (|201ll ). who 
compared the latest theoretical parameter space of the r- 
mode instability against the observed population of accret- 
ing neutron stars in low mass X-ray binaries. 

The present paper makes another contribution to our 
understanding of the /-mode instability and its damp- 
ing. We focus on the action of bulk viscosity, which ac- 
counts for the energy drained from the oscillating fluid as 
a res ult of the induced departure from local beta equilib- 
rium (ISawveJ 19891 : IShapiro fc Teukolskvlll983l ). The novehy 
of our work is the study of bulk viscosity in the nonlinear 
regime (with respect to the oscillation amplitude), thus ex- 
tending previo us models based on amp li tude-indep e ndent 



viscosity (e.g. Ilpser fc LindblomI Il99ll : iLindblomI 1 19951 : 



iGaertig et alj|201ll ) 

The bulk viscosity formalism used in our analysi s is 
based on the work by Alford and collaborators ( Alford ct al.l 
I2OIOI . I201ll . l2012l '). These authors have recently studied the 
role of nonlinear bulk viscosity in the context of the r- 
mode instability (we also briefly discuss the r-mode below 
in Sec I4.3|l . A key result in these papers is an analytical 
expression for the bulk viscosity coefficient in the nonlin- 
ear regime which can be easily incorporated in the standard 
machinery for calculating viscous mode damping. This lat- 
ter procedure consists of solving the inviscid hydrodynamics 
equations and determining the mode frequency and eigen- 
functions, and subsequently using these results in the avail- 
able expressions for the mode's volume-integrated energy 
and damping rate (see section [Sjl. 

Our /-mode calculation is numerica l and is based on 
a 2D time-domain code developed by iPassamonti et al] 
l|2009ah . The code evolves the linearised hydrodynami- 
cal equations of rapidly rotating stars in Newtonian the- 
ory. The frequency and eigenfunctions of specific modes 
can be then extracted from the time-evolved pe r turba - 
tion s with a method developed bv IStergioulas et al.l (|2004| ) 
and iDimmelmeier et al.l l|2006D . An overview and a sample 
test of the code is provided in section [3] In this work we fo- 
cus on the sextupole {I = m = 4) /-mode, which is known to 
be the dominant u nstable multipo l e both in Newtonian and 
Relativistic stars j|Friedman|[l983l : Ilpser fc LindblomI Il99ll : 
iGaertig et aUbOlC " We calculate the resulting instability 
"window" (that is, the unstable portion of the parameter 
space described by the stellar rotational frequency and tem- 
perature) after accounting for the action of shear and bulk 
viscosity (section 14. ip . Our calculation provides results for 
the /-mode saturation amplitude as determined by the bal- 
ance between nonlinear viscous damping and gravitational 
radiation-driven growth (section l4.ip . We also calculate the 
gravitational wave strain associated with an unstable /- 
mode and with its amplitude limited by nonlinear bulk vis- 
cosity (section 14. 2p . Our conclusions and a discussion of re- 
maining issues can be found in section [5] 



eigenfunctions) of the inviscid system. In the second step, 
these results are inserted in the volume integral expres- 
sions for the damping/growth rates. This procedure relies 
on the implicit assumption that the mode evolves quasi- 
adiabatically and the dissipation and growth rates are much 
lower than the oscillation frequency. 

The relevant linearised hydrodynamics equations for 
non-superfluid matter are the familiar Euler, mass conser- 
vation and Poisson equations. In a frame rotating with the 
stellar angular frequency — Qz^ these are: 

dtSv' + 2e''''nj5vk + V (5h + 5$) = -D' , (1) 

P 



dt5p + \/,ip5v') = 0, 



V^5$ = inGSp, 



(2) 



(3) 



where Sv^, Sp and 5$ are, respectively, the perturbed ve- 
locity, mass density and gravitational potential. We have 
also introduced the perturbation of the specific enthalpy, 
5h — SP/p, where SP is the fluid pressure perturbation. The 
term D' in the Euler equation lU is the collective viscous 
force and will be discussed below. 

The system of equations (O-® is closed once an equa- 
tion of state (EoS) is provided. In this work we consider a 
polytropic EoS 



P = Kp'' 



(4) 



where if is a constant and 7 is related to the adiabatic index 
by Af = 1/(7-1). 

Within the present single-fluid hydrodynamics model 
the viscosi ty force consists of shear and bulk viscosity con- 
tributions (|lpser fc LindblomI [T991I ). 



(5) 



where r] and ( are, respectively, the shear and bulk viscosity 
coefficient and 

Sa" = i f\/'5v' + \/'Sv' - ^g''\/5a\ , (6) 



(5(7 = V,5t 



(7) 



where gij is the metric tensor. 

T he viscosity force leads to the following dissipation 
rates l|lpser fc Lindblomlll991^ ■ 

dtEs = 2 j dVT)8(j''5a*j , dtEb = j dVCSaSa* , (8) 

where the integrals are taken over the stellar volume and 
an asterisk denotes a complex conjugate. The gravitational 
radiat ion power is given by the standard multipole for- 
mula (|Thornelll980h . 

dtE^.^ ^ ^Y^Niiuj ~ mnf+^ {\5Di^\' + \5Jim\'') , (9) 



2 FORMALISM 



The strategy for calculating the viscous damping and the 
gravitational wave-driven growth of os cillation modes is 
a we ll established two-step procedure l|lpser fc LindblomI 
[l993) and therefore it will only be outlined here. The first 
step consists of calculating the mode properties (frequency, 



where to is the mode frequency as measured in the rotating 
frame. This expression features mass and current multipoles, 
denoted as SDi,^ and SJim respectively, and the coupling 
constant 



4tvG {I + 1) (/ + 2) 
c2'+M(/-l)[(2/ + l)!!]' 



(10) 
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The /-mode oscillation radiates mainly through the mass 
multipole moments. These are given by 



dVr'SpK 



(11) 



where Yim is the standard spherical harmonic function with 
indices {I, m). 

The above damping/growth rates can be readily calcu- 
lated once the mode solution has been obtained from the 
inviscid Euler equation, i.e. equation ((ij with = 0. The 
inviscid equations also lead to the conserved mode energy 



E : 



dV 



p5v5vi + - {5p5U* + 5p*5U) 



(12) 



where 5U = Sh + 5$. 

With viscosity included, the mode acquires a complex- 
valued frequency. For example, the perturbed velocity has a 
time profile 

iiijt — t/T 



5v^ 



(13) 



and similarly for the other parameters. As a result the mode 
energy evolves according to 



dtE = . 

T 



(14) 



The net combin ed effect of damping an d growth can then 
be expressed as (jlpser fc Lindblomlll99lf ) 



— + — 

Ts -Tb ' 



where 



2E 



dtE, 
2E 



2E 



(15) 



(16) 



are the timescales associated with, respectively, gravita- 
tional wave growth, shear and bulk viscosity damping. A 
mode goes CFS-unstable when r < 0, which requires, as a 
necessary condition, a negative flux dtEf,„ < 0. This can 
happen only when u (cj — mQ) < 0, which is the mathe- 
matical statement of the mode changing from retrograde to 
prograde. 

Up to this point our discussion has been essentially a 
repetition of well known theory. The following sections, how- 
ever, offer a new analysis of the effect of bulk viscosity on 
the /-mode instability. On the other hand there is nothing 
new here about shear viscosity, which anyway plays a minor 
role in our calculation. This is because shear viscosity has 
a significant impact on the instability in the low tempera- 
ture regime (T 10® K) where the instability is anyway 
likely to be s uppressed by superfluid vortex mutual fric- 
tion (|Lindblom & McndcU 1995). Our modelling does not 
account for superfluidity. Without superfluidity the shear 
viscosity is d ominated by neutron collisions and the resulting 
coefficient is l|Flowers fc Itohll 19791 : llpser fc Lindblomlll99ll i 



OAT 9/4rrn-2 -1 -1 

77 = 347p ' T g cm s 



(17) 



We use this coefficient for the shear viscosity in equations ((Sj 
and (fT6l). 



2.1 Nonlinear bulk viscosity 

Considering a neutron star core composed of npe matter, the 
beta equilibrium is established through the modified Urea 



(mUrca) reactions 
N+n ^ N+p+e^ + i?e 



N+p+e ^ iV + n + i^e , (18) 



where A*' is a spectator nucleon that guarantees over- 
all energy-momentum conservation. In more massive neu- 
tron stars where the proton fraction is sufficiently large, 
Xp > 1/9, the more efficient direct Urea (dUrca) process 
(i.e. the re action (|18p without A'^) becomes energetically 
favourable (|Lattimer et al.lll99lh . Although both processes 
are discussed in the following, most attention is given to the 
mUrca case which is the relevant one in a 'canonical' neutron 
star. 

The coefficient of the bulk viscosity associ ated with beta 
reactions takes the form (|Haensel et al.|[200ll ) 



(19) 



where uj is the oscillation frequency in the rotating frame, 
and r represents the particle reaction rate (see Table [l}. The 
temperature power-law is reaction-dependent and takes the 
value g = 3 (g = 2) for mUrca (dUrca). The quantity C is 
the so-called strong interaction susceptibility 



C: 



dpA 



(20) 



where pA = Pn — Pp — pc is the difference between the 
particle chemical potentials on the two sides of ([18} . For the 
simple case of neutron matter modelled as a gas of free (non- 
inter acting) hadrons we have in natural units (|Alford et al.l 
I2OIOI) 



C = 



2/3 



6 m' 



i/3 



(21) 



where mn is the neutron mass. The resulting viscosi ty coef- 
ficient for this type of matter has been calculated bv lSawverl 
(|l98i) and is the most widely used coefficient in the litera- 
ture: 



Co = DxlO p u! Tg cm S 



(22) 



Strictly speaking, the bulk viscosity formulae (|19|) and 
are valid for a 'small' oscillation amplitude in the so-called 
sub-thermal regime, where pA «C fcsT. However, if the os- 
cillation amplitude is sufficiently large the system may en- 
ter the supr a -thermal reg ime pA > ksT. As discussed in 
lAlford et all (|201G| . boill ). this condition may be actually 
reached at relatively small amplitudes Ap/p ^ 0.01. The 
Lagrangian density perturbation Ap is defined in terms of 
the displacement vector as Ap = Sp + ^^Vip. 

In this nonlinear bulk viscosity regime, the formula 
( is replaced by the amplitude-dependent expression 
Alford et al.ll201Gl . l201lh . 



C = Co 



1 + E 



(2i + l)!!x. (CApV^ 
2^(j + l)! \T p ) 



(23) 



where Xj are matter/reaction-dependent parameters given 
in Table [T] This expression is approximate, but remains ac- 
curate as long as the amplitude is not too large. For mUrca 
(dUrca) bulk viscosity th e relevant amplitude thr eshold is 
Ap/p ^ 1 (Ap/p ^ 0.1) (lAlford et al.ll2010l . l201lh . For _an 
oscillation amplitude above these thresholds equation 
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Table 1. This table provides the parameters appearing in the 
expression l|23| l for the the bulk viscosity coe fficient C- Note tha t 
these quantities are given in natural units fsee lAlford et al.l201Cll) . 



Process 


9 


f [MeV^- 


-29] 










mUrca 


3 


4.68 X 10-1^ ( 


' ajpTi ^ 
. "0 , 


J 1/3 






1835^6 
3 ^ 


dUrca 


2 


5.24 X 10-15 ( 


' XpU ^ 

. "a , 


J 1/3 










becomes inaccurate, failing to capture the nonlinear satura- 
tion of the supra-thermal bulk visco sity and the subseq uent 
decreasing behaviour (see fig ure 1 in lAlford et"ai]|201ll ). For 
the purpose of this paper, where the fluid perturbation re- 
mains in the linear regime (i.e. Ap/p < 1), the approxima- 
tion (|23p remains accurate and, in fact, includes the maxi- 
mum viscosity attained in the supra-thermal regime. 

The nonlinear coefficient (|23p provides the main micro- 
physical input for our /-mode instability calculation. A large 
amplitude /-mode is likely to 'activate' the nonlinear terms 
in (I23|l and suffer an enhanced viscous dissipation. If the 
resulting damping rate is sufficiently high then the mode 
may actually saturate and stop growing at some maximum 
amplitude. In the following sections we explore to what ex- 
tent this scenario could be relevant for unstable /-modes of 
rapidly rotating neutron stars. 

In the following, C, is calculated assuming the simple 
case of a free hadron gas, see Table [T] for the relevant pa- 
rameters. For the proton fraction (which is assumed uniform 
throughout the core) we take a typical value = 0.05 when 
mUrca is active and a higher value Xp = 1/9 in the case of 
dUrca viscosity. 

Given the functional form of C, it is also necessary to de- 
fine an appropriate mode-amplitude. We introduce the am- 
plitude parameter q defined by the following expression: 

&vl^ = Q i^R^^Yi^ (7r/2, 0) , (24) 

where R^q is the equatorial radius of the star, and &v^^q is the 
radial component of the velocity perturbation determined at 
r = 7?eq and (Q,<i>) = (7r/2,0). Note that our definition of 
th e mode ampl i tude i s similar to the r-mode definition used 
by I Owen et al.l l| 19981 ). provided that tj is replaced by SI. 

3 NEUTRON STAR MODEL 

Our analysis is based on a Newtonian polytropic stellar 
model with uniform rotation. More specifically, we study 
two polytropes with indices A'' = 1 and = 3/4, which 
provide a reasonable approximation to realistic equations 
of state. The axisymmetric equilibrium configurations are 
determined by solving the stationa ry equat i ons w ith the 
self-consistent method developed by iHachisul l| 19861 ). A de- 
tailed discus sion of t he num erical techniques can be found 
in fPassamonti et al.l (|2009al lbl). In the following calculations 
we fix the stellar mass at M = 1.5 Mq and the constant 
K (in cgs units) at A" = 6.637 x lO" {K = 1.3346) for the 
TV = 1 (A^ = 3/4) model. The resulting stellar radius of the 
nonrotating model is then R — 12.533 km (R — 14.245 km) 



for N — 1 (A' = 3/4 ) . As discussed by llpser fc LindblomI 
(|l99ll ) and iLindblomI l| 19951 ). the /-mode damping/growth 

timescale in neutron stars with different mass and radius 
can be obtained by a simple rescaling. 

Although our /-mode calculation assumes a fluid star 
we implicitly account for the presence of a solid crust by lim- 
iting the action of bulk viscosity in the liquid core, where 
the reactions (jlSp take place. This means that the volume 
integral in equation ((S)) is taken only over the core. This 
truncation should be accurate below the crust melting tem - 
perature r^oit ~ (5-7) x lO'' K (|Shapiro fc T cukolskv IQSi). 
a region which overlaps with the temperature range where 
the /-mode instability is likely to operate. In our model the 
"crust" occupies the outer layer of the star with its boundary 
located at the transition density pec ~ 1.2845 x 10^* g cm"^. 

A step of the calculation that deserves some special dis- 
cussion concerns the constraints imposed on the /-mode am- 
plitude by the chosen stellar models. A first requirement is 
that the mode velocity in the inertial frame should not ex- 
ceed the speed of light. In the rotation range relevant for 
the instability of the I = m — 4 /-mode (which is the most 
unstable multipole), we find that this condition is satisfied 
when a < 1.2 (a < 1.4) for the A^ = 1 (A^ = 3/4) model. 

Another physically motivated (albeit less rigorous) con- 
dition is imposed by the linear perturbation formalism used 
here. We can quantify this requirement in terms of the La- 
grangian density perturbation Ap, as our perturbation for- 
malism (section [21) should break down when Ap/p > 1. This 
condition also limits the use of the approximate bulk vis- 
co sity coefficien t given in equation (|23|) . Indeed, as detailed 
m lAlford et al.l (|201ll ). this approximation for the mUrca 
(dUrca) bulk viscosity is accurate only when Ap/p <. 1 
(Ap/p <. 0.1). In other words, we need to make sure that 
Ap/p remains within the limits of linearised theory. 

The quantity Ap/p is shown in figure [T] for the I — m = 
4 /-mode with a fiducial amplitude a — 1. The density eigen- 
function is calculated at the equatorial plane {6 = 7r/2), and 
in two radial locations, namely, the surface and the crust- 
core interface. All the rapidly rotating models shown in fig- 
ure [T] have Ap/p > 1 at the surface. On the other hand, 
the perturbation remains 'well-behaved' at the crust-core 
boundary, i.e. Ap/p < 1. The same is true for every point 
inside the stellar core. Therefore our analysis can be safely 
used for studying the bulk viscosity in the core. 

Another way of gauging the size of the mode ampli- 
tude is by comparing the mode energy E defined in equa- 
tion H12p against the bulk rotational energy of the star -Erot- 
It is reasonable to expect that in the instability regime we 
have E < E-^ot, as the latter energy is the one tapped by 
the instability. In figure [T] we show these two quantities for 
an /-mode with a = 1. Near the maximum rotation rate, 
the mode energy is much smaller than the rotational en- 
ergy, E ~ (lO"'^ — 10-^)i5rot, but becomes a significant 
fraction of it at relatively lower rotations. For instance, the 
condition E < E^ot is violated in the A^ = 3/4 model when 
f2 < 0.91f2K. However this matters little for our calcula- 
tion since, as we will see in section 14.11 these stars lie out- 
side the instability window of the /-mode. In summary, our 
/-mode calculations preserve self-consistency with respect 
to the mode amplitude (|24|) and the approximate supra- 
thermal bulk viscosity coefficient (|23l) . 
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Figure 1. The energy (left panel) and the amplitude of lS.pl p (right panel) of the Z = m = 4 /-mode in rapidly rotating models with 
C > 0.9 f^K- The stars have a constant mass M = 1.5 Mq and polytropic index A'^ = 1 and A' = 3/4. In all cases, the mode amplitude is 
o = 1. The density perturbation Ap/p is computed on the equatorial plane (Q = 7r/2) at two different radii, at the bottom of the crust 
(r = Rcc) and at the star's surface (r = Req). 



3.1 Numerical Code 

The evolution of the perturbation equations ([l|-(l3| is a 
three-dimensional problem in space. In order to have an 
easier numerical implementation, we can exploit the back- 
ground symmetry and expand the non-axisymmetric per- 
turbations in a Fourier series with respect to the azimuthal 
harmonic ind ex m. As a result, the problem becomes two- 
dimcnsional (jPapaloizou fc PringleH igstf: Pas samonti et al.l 
[2009a, b). For the mass density perturbation the expansion 
is given by 

m — oo 

Sp(t,r,6,(j}) = [5 p'^ (t,r, 9) cos m(l> + 6 (t,r, 9) sinrrKj}] , 

(25) 

and similarly for the other variables. For any m, we evolve 
a system of ten partial difTerential equations for the ten 
variables (5p*, 5$*), while the enthalpy perturbation, 

5h^, is determined by an EoS. For a polytropic star it is 
given by 



on — c, — . 



where is the speed of sound, 



dP 
dp 



= K-yp 



7-1 



(26) 



(27) 



The linearised equations a r e evo lved numerically with a 
code de veloped bv ljones et al.l (|2002l ) and lPassamonti et al.l 
(|2009al lbl). where the reader can find all the details. From 
the time evolution of each perturbation variable, we de- 
termine the mode frequency with a Fast Fourier Trans- 
formation (FFT) and then ext ract the eigenfu n ctions 
with the method deve l oped by IStergioulas et all ()2004l ) 
and IPimmelmeier et al.l ()2006l ). Both the frequencies and 
eigenfunctions are used in section 14.11 to determine the 



damping time of the bulk and shear viscosity and the growth 
time of the gravitational-wave-driven instability. 

3.2 Boundary Conditions 

The time evolution of non-axisymmetric oscillations requires 
the prescription of boundary conditions. We study rotating 
stars with equatorial and rotation axis symmetry, where the 
two-dimensional numerical grid extends over the region < 
r/R{9) < 1 and < 61 < 7r/2. The quantity R{e) is the star's 
radius at a given latitude 9. At each time step, we must then 
specify the variables at the surface, origin, rotational axis 
and equator. 

At the origin (r — 0) and the rotational axis (9 = 0), the 
perturbation variables and equations must be regular. For 
the m > 2 non-axisymmetric modes the regularity condition 
leads to the following equations: 



and Sv^ = , 



(28) 



0, 



(29) 



5p = 5h^5^ = 0, 

both at r = and 9 = 0. 

At the equator {6 = 7r/2), the reflection symmetry di- 
vide s the perturbations into two sets with opposite par- 
ity jPassamonti et cil]|2009al ). For the Z = m = 4 /-mode 
the parity conditions read 

dSp _ dSh _ 95$ 

'W ^ ~d9' ^ ~df 

dSv- e dSv^ 

^ = '^" 

We impose at the surface, r = R{6), the vanishing of 
Lagrangian enthalpy perturbation, 

Ah = 5h + C^^h = 0, (31) 

which is the standard free surface condition used in stellar 
oscillations. The vector field is the Lagrangian displace- 
ment, and the value of the perturbed enthalpy Sh at the 
surface is determined at each time step from equation (|3ip . 
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Figure 2. The left panel shows the I = m = 4 /-mode frequency (in the rotating frame) for our polytropic rotating models. The 
horizontal axis shows the angular velocity in the units of the maximum rotation rate Hk- The right panel compares our results with 
previous work in literature, it displays the viscous and gravitational radiation damping/growth times of the I = m = 4 /-mode for a 
sequence of polytropic rotating stars with N = 1, M = 1.5Mq and T = 1 0^ K . The vertica l axis s hows the various damping times in 
seconds (see legend). Our results (solid lines) and the data determined bv llpser fc LindblomI l ll991^ : [Lindbloinl lll995l ) (symbols) arc in 
good agreement, especially in the CFS instability region of the /-mode. 



3.3 Numerical tests 

W e test our nume r ical f ramew ork using the result s 
of llpser fc LindblomI (|l990l . Il99ll ) and iLindblomI l|l995h . 
These papers study the /-mode instability for various poly- 
tropic Newtonian models, and use an eigenvalue method to 
extract the frequencies and eigenfunctions of the /-mode. 

First of all, we determine the frequency of the onset 
of the CFS instability for inviscid stars. This corresponds 
to models with vanishing /-mode frequency in the inertial 
frame, i.e. u - mO, = 0. We find D. = 0.5604r2o for iV = 1 
and n = 0.5498f2o for A*' = 3/4. These values agre e to bet- 
ter th an 0.5% with the results reported injlpsor & LindblomI 
1I1990I). Note that Op = -v/tt Gpo is the frequency unit used 
bv llpser fc LindblomI (Il990l). where po denotes the average 
mass density of the non-rotating model. In figure [2] (left 
panel) we show the I = m = 4 /-mode frequencies deter- 
mined with our numerical code. 

During the tests we have noticed a slight difference be- 
tween the maximum rotation rat e determined with our nu - 
merical code and that reported in llpser fc LindblomI l|l99lf ). 
Although the difference is very small (less than 0.6%) it 
has a visible impact on the instability window. For the 
iV = 1 and iV = 3/4 polytropes, we fin d = 0.6352^0 
and O k = 0.6454S7o, respectively, while llps er fc LindblomI 
l|l99ll) report D,k = 0.639r2o and Qk = 0.64800- We have 
no reason to question the reliability of our code, given that it 
reproduces to hi gh accuracy th e axisymmetric background 
configurations of iHachisul (|l986l ). Hence, we think that this 
small discrepanc y is due to the differ ent numerical tech- 
niques used in Ipser fc LindblomI (|l99ll ). 

As a final test, we compare the damping/growth times 
of the gravitational radiation a nd the viscous damping 
times with the data extracted from llpser fc LindblomI (|l99ll ) 
and lLindbloinl (|l995l ). The results are shown in figure [2] for 



the A*' = 1 rotating models with te mperature T = 10 K. We 
have rescaled the data reported in llpser fc LindblomI (|l99ll ) 
with the flK of our equilibrium configuration. By using the 
standard bulk viscosity given in equation (|22|l . we find a 
very good agreement with previous results. 



4 RESULTS 

We are now ready to present our results on the impact of 
nonlinear bulk viscosity on the gravitational wave-driven /- 
mode instability. We first discuss how the mode's instabil- 
ity window is modified as a result of viscous damping (sec- 
tion I4.1|l . We then calculate the gravitational wave strain 
associated with the /-mode and provide upper limits for 
its detectability by present day and future detectors (sec- 
tion I4.2|l . In the final part of our calculation (section I4.3|l 
we make a brief digression and discuss the role of nonlinear 
bulk viscosity in the context of the r-mode instability. 



4.1 Instability window 

The boundary of the /-mode instability window on the T 
plane is determined by 

i = . (32) 

T 

For a given T this condition provides a critical rotational 
frequency ilc above which the mode is CFS-unstable. 

The previous studies of the /-mode instability included 
bulk viscosity in the su b-thermal reg ime, typically using 
Sawyer's coefficient (|22[) (|Sawveilll989l ). In that case equa- 
tion (|32|l . and consequently Qc, are independent of the mode 
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Figure 3. This figure displays the instability window of the I = m = A /-mode for the N = 1 (left panel) and TV = 3/4 (right panel) 
polytropic models with mass M = 1.5 Mq. The bulk viscosity is determined using the coefficient f (equation 1231 1) for the modified 
Urea process. For various amplitudes a we show the effects of the nonlinear bulk viscosity on the instability window. The shaded regions 
represent the temperature range where neutron superfluidity (and proton superconductivity) in the core is expected to be present. The 
critical temperature Ten ~ 5 X 10^ K for the onset of neutron superfl udity in the core represents a lower l imit, as suggested by the 
recent models for the observed cooling of the Cassiopeia A neutron star ijPage et al.ll201ll : IShternin et al.|[2011^ . Below that temperature 
superfluid mutual friction is expected to suppress the instability. 



amplitude a. However, once the full nonlinear bulk viscos- 
ity expression ((23} is used, the instability window depends 
explicitly on a. 

We first consider bulk viscosity due to the mUrca pro- 
cess. In figure Owe show the / = m = 4 /-mode instability 
window for our A*' = 1 and N = 3/4 polytropic models and 
for different mode amplitudes. The effect of nonlinear bulk 
viscosity is clearly more pronounced in the N = 1 model, 
and for amplitudes a > 0.1. For an amplitude as large as 
Q ~ 1 the window area is decreased to a significant degree, 
resulting in a minimum Sic ~ 0.975. The stiffer TV = 3/4 neu- 
tron star model is less sensitive to nonlinear viscous damp- 
ing. This is not surprising given that the instability's growth 
time Tgw is shorter than the one in the A'' = 1 polytrope. To 
summarise, we can conclude that nonlinear bulk viscosity 
due to mUrca reactions has a moderate effect on the /-mode- 
instability even for the largest mode amplitudes consistent 
within our framework. 

The instability window in models with dUrca reactions 
is shown in figure|4] As expected, the higher efficiency of the 
dUrca reactions leads to a stronger bulk viscosity which sup- 
presses the instability above T « 10^ K. This temperature 
is very close to the expected onset of neutron superfluid- 
ity in the core, and once neutron superfluidity is present 
the inst ability is suppressed as a result of vortex mutual 
friction ^Lindblom fc MendellllTgQSl ). The resulting /-mode 
instability window is much smaller than the (already small) 
window of the mUrca models. Focusing on the damping of 
nonlinear bulk-viscosity we find that the situation is similar 
to the previous mUrca models. The critical frequency fie is 
moderately increased, more for the A' = 1 polytrope and 
less for the A'' = 3/4 one. 

Our results suggest that bulk viscosity, despite its 



increased efficiency due to the nonlinear supra-thermal 
contribution, can only cause a moderate degree of ad- 
ditional damping. How robust is this conclusion? Based 
on re cent results on the /-mode instability in relativistic 
stars (jGaertig et al.|[201ll ) which generally predict a shorter 
growth timescale than the Newtonian one ~ O.lr^""*), 
one would suspect that the role of nonlinear bulk viscosity 
is even less important. At the same time, however, we also 
need to account for the fact that bulk viscosity itself be- 
comes stronger, roughly by a factor ~ 10, as we move to a 
more realis tic model with interacting hadrons (jAlford et al.l 
12011 12OI2I ). With these two effects combined and the en- 
suing balance between them (at least approximately) it is 
likely that our Newtonian analysis leads, after all, to robust 
results. Obviously, a definitive answer can only be reached 
by means of a more detailed calculation (e.g relativistic stel- 
lar models with a more realistic equation of state). 

As we have already pointed out, our analysis ne- 
glects the presence of neutron and proton superfiuidity. 
At the same time, however, we know that superfluidity 
plays a key role by suppressing the /-mode instability in 
sufficiently cold neutron stars through vortex mutual fric- 
tion (Lindblom & McndcU 1995). Bulk viscosity is also di- 
rectly affected by superfluidity as a consequence of the re- 
duced rates of the mUrca/dUrca reactions in superfluid neu- 
tron stars. Thus we have good reasons for revising the bulk 
viscosity damping in superfluid neutron stars. Within our 
framework we can do that by introducing the appropriate 
superfluid correction factor in the viscosity coefficient (see 
below). Other important aspects of superfluid hydrodynam- 
ics, like damping due to vortex mutual friction and other vis- 
cous degrees of freedom (|Andersson fc Comerll2006h . can not 
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Figure 4. This figure siiows ttie instability window of the Z = m = 4 /-mode for the N = 1 (left panel) and TV = 3/4 (right panel) 
polytropic models with mass M = 1.5 Mq. Differently from figureO the bulk viscosity is here generated by the direct Urea process. See 
figure [3] for more details about the notation. 



be accessed wi thout the full machinery of multi -fluid hydro- 
dynamics (e.g. lAndersson. Glampedakis fc Has kell 2009). 

As a newborn neutron star cools down the transi- 
tion to proton superconductivity is likely to occur first, 
with neutrons in the core becoming superfluid at a later 
stage. This orde r is the one sug gested by pairing calcu- 
lations (e.g. Ka minker et al.ll2002l ). The recent models for 
the observed thermal evolution of the Cassiopeia A neu- 
tron star provide constraints for the critical temperatures 
for the onset of proton and neutron superfluidity in the core, 
Tcp ^ (2 - 5) X 10^ K and T en ~ (5 - 9) x 10* K (|Page et all 
I2OIII : IShternin et al] I2OIII ). There is therefore a tempera- 
ture range where the stellar core contains normal neutrons 
and superconducting protons. In this regime, where vortex 
mutual friction cannot yet operate, it makes sense to study 
how bulk viscosity is modifled by proton superconductivity. 

This can be achieved by using the modifled viscosity 
coefficient 

C=up = C7^, (33) 

where is the bulk viscosity coefficient for normal mat- 
ter defined by equation (|23p . The factor TZ represents the 
superfiuidity-induced reaction rate suppression and is a 
function of the critical temperature Tcp for the onset of pro- 
ton supercondu ctivity. The explicit form of TZ is given by 
equation (34) in lHaensel et al.l (|200ll ). It should be pointed 
out that the modification (|33p is the appropriate one for the 
bulk viscosity coefficient in the sub-thermal regime. As far 
as we know there is no similar result available for the supra- 
thermal hence the superfiuid modification H33|) should be 
viewed as an approximation. 

The impact of proton superconductivity on the /-mode 
window is shown in figure [5] We see that proton supercon- 
ductivity can mitigate to a significant degree the damping 
action of nonlinear bulk viscosity. This effect becomes more 
prominent with an increasing critical temperature Tcp. 

Nonlinear bulk viscosity, when efficient, is a mech- 



anism for saturating a growing mode due to its intrin- 
sic dependence on the mode-amplitude. Of course this 
does not exclude the simultaneous presence of other com- 
peting (and perhaps more efficient) saturation mecha- 
nisms. A strong candidate alternative saturation mecha- 
nism is nonlinear mode coupling. For instance, mode cou- 
plings with other inertial modes is known to saturate the 
r-mode instability at small amplitudes , irrespectively of 
rotat ion/temperature (|Arras et al ] l2003l : iBondarescu" et al.l 
I2OO7I ). Thus it is also possible that a similar effect could 
be at work in the case of the /-mode. Up to now 
there is limited work on the nonlinear hydrodynamics 
of the /-mode (lOu et al.1 12004 IShibata fc Karinol |2004| : 
iKastaun et al.ll2010l ). and none for the / = m = 4 multi- 
pole which is of interest to us. 



Still, we can obtain some hints from the available work 
on the quadrupole /-m o de. T he nonlinear numerical simula- 
tions of iKastaun et all (|201Cll ) suggest that the /-mode am- 
plitude in relativistic polytropic stars is genetically limited 
by wave breaking at the surface and mode couplings. Con- 
sidering a rotating A'^ = 1 polytropic model with axis ratio 
Rp/R eq = 0.7, where Rp is the polar radius, IKastaun et al.l 
(|2010l ) found that an Z = 2 /-mode with an initial energy 
i? ~ 2 X 1O~^M0C^ is efficiently damped by nonlinear ef- 
fects on a timescale of about 10 ms. Assuming that this 
result remains valid also for the m — i /-mode we can de- 
duce (with the help of figure [T] and noting that the above 
axis ratio corresponds to ~ 0.9Qk) that nonlinear mode 
couplings could saturate the mode amplitude at a < 1. This 
exercise suggests that nonlinear hydrodynamics may provide 
the dominant saturation mechanism for the /-mode insta- 
bility. However, without a detailed analysis of the I = m = 4 
/-mode dynamics we cannot yet draw any secure conclusion. 
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Figure 5. This figure shows the impact of proton superconductivity on the instability window of the I = m = A /-mode for different 
choices of proton critical temperature Tcp (assumed to be uniform throughout the core). The stellar mass is M = 1.5 Mq and the 
polytropic index N = 1. The mode amplitude is a = 0.1 (dashed lines) and a = 0.5 (solid lines) in the left panel, and a = 1 in the right 
panel. We show in the right panel also the critical curve for the sub-thermal bulk viscosity (dashed-line) . 



4.2 Gravitational wave signal 

The final part of our /-mode instability analysis is devoted 
to the calculation of the gravitational wave st rain and its 
detectability. Unlike previous studies IjLai fc Sh apiro 1995; 
lOu et al.|[200i : IShibata fc Karinol |2004| ) we concentrate on 
the I = m = A /-mode which is, as we have already pointed 
out, the most instabili ty-prone multipole in both Newtonian 
and relativistic stars dCaertig et al.ll2011^ . 

The gravitational wave strain can be computed by 
means of the standard multipole formula (|Thornelll980l ) 



44 _ G 1 (fx 



(34) 



where I*'' is proportional to the mass moment (|lip . 

2^"" = It7 ^'^■044, (35) 
945 

and Tij' is the corresponding s pin tensor h armonic angular 
function of the "electric- type" (|Thornelll98Q, ). 

It turns out that the two independent wave polariza- 
tions can be expressed in terms of the spin-weighted spher- 
ical harmonic -2Y 



7,44 -7 44 7.44 -(^44 

hge - ihg^ = h -2Y , 
where the wave amplitude ft^* is defined as 
/To 471 G 



945 r c6 dt4 

For a monochromatic wave equation p7l) reduces to 
, 44 47r G 



945 r c' 



4 



dV5p{r) rYl^e' 



he 



(36) 



(37) 



(38) 



where uji = u — mO. is the mode frequency in the inertial 
frame. 

As a representative gravitational wave amplitude wc nu- 
merically calculate h'^*' using the relevant /-mode frequency 



and eigenfunctions. In figure [5] (left panel) we show our re- 
sults for h assuming a source located at a distance of 20 Mpc 
and a mode amplitude a = 1. As expected, the gravitational- 
wave signal is stronger for models with Q. ~ Q,k, where the 
star is more unstable and uji is higher. 

For assessing the actual detectability of a neutron star 
undergoing oscillations due to an unstable /-mode we need 
to take into account the cumulative effect of the observed 
number of oscillation cycles in a given frequency bandwidth. 
This information is encoded in the characteristic wave strain 



(39) 



where v = u)i/2Ti and (...) denotes the average over the 
angles {9, (j)). 

The time profile of the mode frequency can be de- 
ter mined following the two-stage evolution of the instabil- 
ity (|Owen et al.l ll9S?^ ). The onset of the instability leads 
to the first (linear) phase of exponential growth where 
u « constant. Eventually this growth arrives to a halt as 
the mode is saturated by nonlinear physics (for instance, 
the nonlinear bulk viscosity studied in this paper). During 
this second (nonlinear) phase the mode amplitude remains 
constant and the emitted gravitational radiation removes 
angular momentum from the star's bulk rotation causing its 
spin-down. This evolution leads to a monotonically decreas- 
ing mode frequency uit). Here we consider only the second 
nonlinear saturation phase and neglect the thermal evolu- 
tion of the star and other mechanisms which can accelerate 
the spin-down. Our results provide therefore an approximate 
upper limit for the gravitational wave strain. A more de- 
tailed analysis of the detectability of unstable /-modes in 
neutron stars will be presented in future work. 

In order to estimate the number of oscillations A'cyc = 
\dv/dt\~^ near the given mode frequency u we need to 
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Figure 6. Gravitational-wave strain generated by the unstable Z = m = 4 /-mode for rotating polytropic models with indices TV = 1 
and N = 3/4. The stars have mass M = 1.5 Mq and are located at 20 Mpc. On the horizontal axis we show the inertial /-mode 
frequency, which can be directly related to the star's rotation rate (see the left panel of figure [2]l. The left panel displays the maximum 
strain h for an /-mode with fiducial amplitude a = 1. In the right panel, we show the gravitational- wave signal emitted by the /-mode 
during the nonlinear saturation phase. We compare the characteristic strain he for different mode amplitudes with the hrms of the 
gravitational- wave detectors (see legend). The two curves with larger he correspond to an ideal case where the detectors can follow the 
entire evolution of the /-mode instability, and he is then amplitude independent. The other curves represent a signal integrated for 1 
week. Note that we have neglected the thermal evolution in our spin-down model, the curves in the right panel may therefore represent 
the gravitational radiation emitted by a star which slows down at constant temperature from the Kepler frequency to the minimum 
of the instability window. For the mUrca process, we specify also the characteristic strain emitted by a model at the minimum of the 
critical curve. For an /-mode with a = 1 (a = 0.1) this signal is denoted by filled circles (triangles). 



(40) 



consider the angular momentum balance for the system, 
dJ _ _ dJ 

dt ^ ~ dt 

where J = IQ, is the steUar angular momentum B Writ- 
ing dJ/dt = {dJ/dfl){dn/dt), dQ/dt = {dQ./ dv){dv / dt) and 
solving with respect to the derivative of the mode frequency, 
we obtain 

dv 1 diu dj f dJ^ ^ 

H ^ 2ndn lit \dn 

The radiated angular momentum is given by l|Thornelll980l ) 

dJi 



(41) 



dt 



= 2N4UJ1 l(5D44i <5i; 



(42) 



where Siz is the Kronecker delta. From this expression it is 
obvious the scaling dJ /dt\g^ ~ , which means that A^cyc ~ 

. Then from (|39|) we can deduce that the characteristic 
strain he is independent of the mode amplitude a. 

The above calculation is idealised in the sense that it 
implicitly assumes a gravitational wave observation lasting 
as long as the /-mod e instability evolu tion timescale. This 
timescale is given by l|Owen et al.lll99^ ) 

dt 



'?~cvol — 



(43) 



^ To be more precise we should also consider the variation of the 
mode's canonical angular momentum, but its contrib ution is gen- 
erally smaller than the stellar angular momentum llOwen et al.l 
Il998h and hence can be neglected for simplicity. 



For a fiducial mode amplitude a = 1 and a. N = 1 
model, we find that for ~ Q.k the evolution timescale 
is Tcvoi ~ 10 weeks, and increases to TcvoI ~ 2 yr for slower 
rotating models with Q. ~ Q.Q^^Ik- The evolution timescale 
becomes even longer for modes with a < 1 given the scal- 
ing Tcvoi ~ oT^ . Thus in the most relevant cases the mode 
evolution could be much longer than the typical observation 
time of LIGO and Virgo. What that means is that from 
the detectors' point of view the signal can be considered as 
essentially monochromatic. Besides the temporal limits set 
by the detectors' capability, there might be additional ones 
from the physical system itself. For instance, the star may 
efficiently spin-down as a result of a strong magnetic field, 
thus quenching the /-mode instability and the emission of 
gravitational radiation. 

Therefore, in order to provide a 'realistic' estimate for 
the gravitational wave strain he associated with the /-mode 
instability we assume a monochromatic signal of duration 
Iweek, for a source at a 20 Mpc distance. For the purpose 
of comparison we also consider the ideal case of a signal 
lasting a time period TcvoI- The results are shown in figure [6] 
together with the sensitivity curves of the current and fu- 
ture generation Earth-based gravitational wave detectors. 
The detectors' noise curves are computed in the standard 
way using /irms = \/ vSh {v), where Sn {f) is the detector's 
power spectral density. In our simplified spin-down model we 
have neglected the effects of the neutron star's cooling. The 
curves shown in figure [6] (right panel) may represent there- 
fore the signal emitted by an unstable star which moves at 
constant temperature from the Kepler frequency down to the 



The f-mode instability in neutron stars 11 



minimum of the instability window (see figure [3|. However 
the thermal evolution may significantly affect the trajectory 
through the instability window and the star may therefore 
leave the instability region at larger rotation rates and thus 
at higher mode- frequencies than those indicated in figure [6l 
Our results suggest that the /-mode instability may be 
detectable from the future ET for a mode amplitude a > 
0.1. Detection by Advanced LIGO/Virgo requires a larger 
amplitude a ~ 1 and stars with a stiff equation of state (in 
the figure the = 3/4 polytrope is clearly favoured over 
the A = 1 one). 



4.3 Digression: the r-mode instability 

So far we have discussed the role of (nonlinear) bulk viscos- 
ity in the damping of the /-mode instability. However, as 
mentioned in the introduction, the physics of nonlinear bulk 
viscosity has alre ady been discussed in the conte xt of the r- 
mode instability (|Alford et al.ll2010l . l201ll . 120121 '). The main 
conclusion of that recent work is that bulk viscosity operat- 
ing in the supra-thermal regime (see section E7I|) could sup- 
press the r-mode instability (regardless of the stellar tem- 
perature), and prevent the growth of the mode amplitude 
above some maximum saturation level. In this section we 
revisit the topic of r-mode nonlinear viscous damping with 
the aim of providing a better understanding of the recent 
results. Why this is necessary will become obvious in the 
following discussion. 

To begin with, we introduce th e 'standard' dimen sion- 
less r-mode amplitude a defined as (|Owen et al.lll998l ) 



rBj iujt 



(44) 



where is the magnetic vector harmonic. Obviously a 

must not be confused with the /-mode amplitude defined in 
equation (|24p . This expression describes the r-mode velocity 
at leading order with respect to the stellar spin in the famil- 
iar sl ow-rotation approximation (see lAndersson fc Kokkotad 
I2OOII . for further discussion). In this approximation the star 
maintains its spherical shape and the mode's energy is dom- 
inated by the kinetic part, i.e. 
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where the I = m r-mode is considered. The calculation of the 
bulk viscosity damping timescale requires the input of the 
density fiuctuation associated with the perturba tion (|44|) . 
After adding a missing factor inlOwen et al. I (|l998t ) (see also 
a comment in iLindblom et al.lll999[ ) 
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where 5$o(r) is the radial eigenfunction of the perturbed 
gravitational potential. At this point it is important to 
make a c larification. I n their study of the r-mode viscous 
damping lAlford et all l|201ll , [ 2012) have used a different 
definition for the r-mode amplitude. More specifically, Al- 
ford et al. use the amplitude a defined in equation (4.1) 



of iLindblom et all (|l999l )FI The relation between the two 
different de finitions of the r-mod e amplitude is (see also the 
footnote in ILindblom eraLlll999l ) 



a ^ a,j — (m + lf (2m + 1)! 
m 



(47) 



For the most uns table m = 2 r-mode we fi nd a ~ 713. 

The fact that I Afford et all (|201ll , l2012h use the notation 
a for the non-standard r-mode amplitude a may be the 
cau se for some confusio n. According to their results, e.g. fig. 
2 m lAlford et al.l (|201ll ). bulk viscosity in the supra-thermal 
regime is capable of suppressing the r-mode instability for 
a wide range of stellar temperatures provided a > 0.05. 
However, in terms of the standard r-mode amplitude (|44|l 
this result corresponds to a saturation amplitude a ^ 3.5. 
Obviously this is a rather large amplitud^l, much higher 
than the saturation amplitude a <. 10~^ due to t he r-mode's 
nonlinear coupling with othe r inertial modes l|Arras et al.l 
I2OO3I : iBondarescu et al.ll2007l ) . Thus, we are led to conclude 
that for an amplitude a < 1 the r-mode instability is affected 
very little by the nonlinear contribution of bulk viscosity. 



5 CONCLUDING REMARKS 

We have studied the gravitational wave-driven instability of 
the /-mode in Newtonian polytropic neutron star models. 
The novelty of our work is the inclusion of bulk viscosity in 
the supra-thermal regime. We have found that the ensuing 
nonlinear enhancement of bulk viscosity leads to a moderate 
change in the /-mode instability window. For the particular 
case of the Z = m = 4 /-mode (the most unstable multipole) 
we have shown that nonlinear bulk viscosity associated with 
the modified Urea process has a notable effect provided the 
mode amplitude, as defined in equation (|24p. is a > 0.1 
(see figure |3} . This amplitude translates to a mode energy 
E > 10~^i?rot for Q « f^K, increasing to i? > 10"^i5rot 
for f2 « f2c. However, even for the extreme case of a mode 
amplitude as large as a ~ 1 bulk viscosity cannot entirely 
suppress the instability. 

If the more powerful direct Urea reactions can oper- 
ate (which would be the case in the most massive neutron 
stars) then the resulting /-mode instability window is sig- 
nificantly reduced already in the sub-thermal bulk viscos- 
ity regime. The transition to the supra-thermal regime has 
again a moderate effect to the shape of the instability win- 
dow. It turns out that in the temperature range where the 
change is important the instability is already likely to be 
suppressed by superfiuid mutual friction (see figure |4)|. 

The presence of superfiuid components in neutron star 



^ Shortly after the submission of this paper, Alford et al. have 
revised their r-mode calculation and used the more suitable am- 
plitude defined by equation 1144 I I. 

This can be demonstrated if we take a. N = 1 polytropic star 
with mass M = 1.4M0 and radius R = 12.533 km, rotating at the 
Kepler limit (the maximal rotation rate can be approximated with 
Ok = ^V'^'dpo, where po is the average density at f2 = 0). We 
then obtain Bkin = 1.126 X lO^^^ erg and Erot = 1 .8 X 10^^ erg 
Hence, the r-mode saturation amplitude obtained in llAlford et al.l 
[2011, 2012) would imply a mode energy comparable to, or even 
larger, than the stellar rotational energy. 
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matter could diminish the impact of bulk viscosity by sup- 
pressing the rates of the chemical reactions involved. By con- 
sidering a temperature regime where the protons are super- 
fluid and the neutrons normal and by treating in an approx- 
imate manner the superfluid suppression, we have shown 
that nonlinear bulk viscosity is significantly weakened (see 
figure [5} . 

Apart from calculating the bulk viscosity-modified /- 
mode instability window we have also provided upper lim- 
its for the gravitational wave signal from a mode with 
a = 0.1 — 1, i.e. the amplitude range where nonlinear bulk 
viscosity can partially saturate the mode (figure [6]). Our re- 
sults suggest that the Z = m = 4 /-mode instability (oper- 
ating in a neutron star at the fiducial distance of 20 Mpc) 
could be detectable by a third generation instrument like 
the ET for an amplitude a > 0.1 and an observation time of 
~ 1 week. Detectability by Advanced LIGO/Virgo requires 
a rather high amplitude (a ~ 1) and/or a longer obser- 
vation time (e.g. the evolution timescale (|43p ). Moreover, 
the /-mode detectability is moderately sensitive to the stel- 
lar polytropic index, with lower indices (stiffer equation of 
state) being favoured. It should be pointed out that our work 
provides only 'static' upper limits, in the sense that the com- 
bined rotational/temperature evolution of an /-mode unsta- 
ble neutron star has not been properly incorporated. A more 
detailed analysis of this issue will the subject of future work. 

Our brief excursion to the study of the damping of the 
r-mode instability by nonlinear bulk viscosity has revealed 
that a large mode amplitude (a > 1, see equation (|47p ) is 
required for the effect to be relevant. Thus, most likely, the 
saturation of the r-mode is controlled by nonlinear mode 
couplings rather than bulk viscosity. A similar statement 
could also be true for the /-mode but this remains to be 
(dis)proved, given the lack of any concrete results on the 
/-mode nonlinear couplings. Some light on this issue can 
be shed by a r ecent study of the n onlinear dynamics of the 
I = 2 /-mode (iKastaun et al.ir2010l) . Extrapolation of those 
results to the I = 4 mode suggests that nonlinear hydro- 
dynamics could outperform bulk viscosity in saturating a 
growing mode. 

Besides the obvious need for a better understanding 
of the /-mode nonlinear dynamics, several other improve- 
ments are required in this research area. For instance, our 
own analysis can be clearly improved by moving from New- 
tonian to General Relativistic gravity and from polytropic 
matter to more realistic equations of state. Indeed, as sug- 
gested by recent work on the relativistic /-mode, the insta- 
bility's growth time may shrink as much as a factor ~ 10 
compared to the Newtonian result. However, as we pointed 
out in section IT2] this effect is likely to be counter-balanced 
by a similar increase in the bulk viscosity strength as a re- 
sult of using more realistic equations of state. The presence 
of a solid crust adds one more level of complexity to the 
problem. Despite its central role in lim iting the r-mode in- 
stability l|Andersson fc KokkotadboOll ). the dynamics of the 
crust has been totally ignored in the context of the /-mode 
instability. Another interesting possibility arises when ex- 
otic matter (for example, deconfined quarks) is present in 
the stellar interior. Previous work on the r-mode instabil- 
ity in neutron stars with exotic cores has unveiled several 
important diff erences with respect to normal hadro nic neu- 
tron stars (e.g. lMadsenlll992l : lAndersson et al.ll2002l ). This is 



likely to be the case also for the /-mode instability. Clearly, 
all the above issues are of great interest and we plan to ad- 
dress some of them in the near future. 
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